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Purpose We present a new morphometric measure of trabecular bone micro-architecture, called mean node 
strength (NdStr), which is part of a newly-developed approach called long range node-strut analysis. Our 
general aim is to describe and quantify the apparent “lattice-like” micro-architecture of the trabecular bone 
network. 

Methods Similar in some ways to the topological node-strut analysis introduced by Garrahan et ah, our 
method is distinguished by an emphasis on long-range trabecular connectivity. Thus, while the topological 
classification of a pixel (after skeletonisation) as a node, strut, or terminus, can be determined from the 3x3 
neighbourhood of that pixel, our method, which does not involve skeletonisation, takes into account a much 
larger neighbourhood. In addition, rather than giving a discrete classification of each pixel as a node, strut, 
or terminus, our method produces a continuous variable, node strength. The node strength is averaged over 
a region of interest to produce the mean node strength (NdStr) of the region. 

Results We have applied our long range node-strut analysis to a set of 26 high-resolution peripheral 
quantitative computed tomography (pQCT) axial images of human proximal tibiae acquired 17 mm below 
the tibial plateau. We found that NdStr has a strong positive correlation with volumetric trabecular bone 
mineral density (BMD). After an exponential transformation, we obtain a Pearson’s correlation coefficient of 
r = 0.97. Qualitative comparison of images with similar BMD but with very different NdStr values suggests 
that the latter measure has successfully quantified the prevalence of the “lattice-like” micro-architecture 
apparent in the image. 

Moreover, we found a strong correlation (r = 0.62), between NdStr and the conventional node-terminus 
ratio (Nd/Tm) of Garrahan et al. The Nd/Tm ratios were computed using traditional histomorphometry 
performed on bone biopsies obtained at the same location as the pQGT scans. 

Conclusions The newly introduced morphometric measure allows a quantitative assessment of the long- 
range connectivity of trabecular bone. One advantage of this method is that it is based on pQGT images that 
can be obtained noninvasively from patients, i.e. without having to obtain a bone biopsy from the patient. 
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I. INTRODUCTION 

Several studies have shown that measures of trabec¬ 
ular bone micro-architecture and bone strength are 
correlated.fl^ Together with loss of bone mass, changes 
in the trabecular bone micro-architecture occur during 
ageing, during development of osteopenia and osteoporo¬ 
sis as well as in connection with immobilisation or space 
flight, and can lead to an increased risk of bone fracture. 
The vertebral bodies and the epiphyses and metaphyses 
of the long bones consist mainly of trabecular bone sur¬ 
rounded by a thin cortical shell.A dramatic change 
in the state of the trabecular bone leads to an increased 


fracture risk.^^ 

Bone mineral density (BMD) is the most commonly used 
predictor of bone strength and fracture risk, and also 
the most commonly used general descriptor of the state 
of the bone. Non-linear relationships have been estab¬ 
lished between volumetric B MD and compressive bone 
strength and elastic modulus.^^^ However, for trabecu¬ 
lar bone, it has been established that a part of the vari¬ 
ation in the strength of the bone cannot be explained by 
BMD alone, but is instead due to the micro-architecture 
of the trabecular network. For example, a relationship 
has been established between the mechanical properties 
of the bone and the shape, orientation, bone trabecular 
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volume fraction, and thickness of the trabecnlaej ^ * ^ * ^^^ ^^ 
A series of new methodologies based on techniques from 
nonlinear data analysis has also been introduced in order 
to study the relationship between the complexity of the 
trabecular bone network and bone strengthP®^ Saparin 
et al. established such a relationship by use of struc tural 
measures of complexity based on symbolic encodingPSEll 
Furthermore, several studies using numerical modelling 
of the trabecular bone and finite element analysis have 
also confirmed the importance of the trabecular bone 
micro-architecture to bone strengthP^H^ 


In the present study we propose a new method for anal¬ 
ysis of trabecular bone micro-architecture from high- 
resolution quantitative computed tomography (QCT) 
images, as at present it is still not possible to perform 
micro CT (/iCT) on humans in vivo and in situ. In con¬ 
trast to classical histomorphometry, CT images can be 
obtained in a nondestructive and noninvasive manner, 
which is preferable in a clinical setting. Our method 
uses a new approach called long range node-strut analy¬ 
sis that quantifies the apparent nodes and struts in the 
trabecular network. In contrast to the node-strut anal¬ 
ysis of Garrahan et al.,1^ our method emphasises long 
range connectivity of the trabecular network, over a dis¬ 
tance controlled by a parameter in the algorithm. The 
analysis has the dual aims of describing the shape of the 
trabecular network and predicting bone strength. 


The trabecular bone compartment consists basically of 
two components: bone and marrow. Due to the limited 
resolution of present day CT and peripheral quantita¬ 
tive CT (pQCT) scanners it is not possible to completely 
resolve the trabecular bone micro-architecture. This re¬ 
sults in variations of the CT values of the trabecular vox¬ 
els, even if the intrinsic bone density is constant through¬ 
out the trabecular network. Consequently, our method 
takes this apparent variation in bone density into consid¬ 
eration rather than segment or binarise the image, which 
would imply a loss of this additional information. We 
assume that higher CT values of a given trabecula will, 
all else being equal, result in a higher compressive bone 
strength. We note that our method does not require 
skeletonisation of the image. 


We apply our method to 2-dimensional pQCT images 
of human proximal tibiae and quantify the trabecular 
bone micro-architecture at different levels of bone in¬ 
tegrity ranging from normal healthy bone to osteoporotic 
bone (as assessed by their BMD). We compare our re¬ 
sults with the node/terminus ratio (Nd/Tm),l^ com¬ 
puted using traditional histomorphometry performed on 
bone biopsies obtained at the same location as the pQCT 
scans. 


II. MATERIALS 

The study population comprised 18 women aged 75-98 
years and 8 men aged 57-88 years. At autopsy, the tibial 
bone specimens were placed in formalin for fixation. For 
each specimen a pQCT image and a bone biopsy were 
obtained from the same location. 

For each proximal tibia, an axial QCT slice was acquired 
17 mm below the tibial plateau with a Stratec XCT-2000 
pQCT scanner (Stratec GmbH, Pforzheim, Germany), 
with an in-plane pixel size of 200 fim x 200 /im and a 
slice thickness of 1 mm. In some cases, the scans were 
performed after the biopsies were taken. Therefore, the 
holes left from the bone biopsy appear in some of the 
pQGT images. A standardised image pre-processing pro¬ 
cedure was app lied to exclude the cortical shell from the 
analysisOne of the resulting images is shown in 


Fig. [T1 



FIG. 1. An axial pQCT slice of human proximal tibia ac¬ 
quired 17 mm below the tibial plateau. The cortical shell 
has been removed from the image. The trabecular BMD of 
this sample is 107 mg/cm^. The hole visible on the left of the 
image is the result of a cylindrical biopsy. 

Cylindrical bone samples with a diameter of 7 mm 
were obtained 17 mm distal from the centre of the me¬ 
dial facet of the superior articular surface by drilling 
with a compressed-air-driven drill with a diamond-tipped 
trephine at either the right or the left proximal tibia. 
These bone biopsies were embedded undecalcified in 
methyl methacrylate, cut in 10-/im-thick sections on a 
Jung model K microtome (R. Jung GmbH, Heidelberg, 
Germany), and stained with aniline blue (modified Mas¬ 
son trichrome). The mounted sections were placed in an 
fiat-bed image scanner and 2540 dpi digital 1 bit images 
of the sections were obtained as previously described in 
detail.l^ The resulting pixel size is 10 fim x 10 /im. 

The trabecular BMD of each pQCT slice was calculated 
using a linear relationship derived on the basis of experi¬ 
mental calibration with the European Forearm Phantom, 
as described by Saparin et al.^^ The trabecular BMDs of 
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the slices in this study range from 30 to 150mg/cm^, 
with a median of 97 mg/cm^. 


III. ANALYTICAL METHOD 


Most of this section describes the new image analysis 
method, long range node-strut analysis^ which includes 
the new measure, mean node strength (NdStr). This 
method was applied to the pQCT sections described in 
Section [Ilj after removal of the cortical shell. We also 
computed a standard measures from the same regions of 
interest of the same pQCT images: trabecular volumet¬ 
ric bone mineral density (BMD), calculated as described 
by Saparin et al.^^ 

For further comparison, topological 2-dimensional node¬ 
strut analysi^^ was performed on the histological sec¬ 
tions described in Section [TT| using a custom-made com¬ 
puter program.l^ The trabecular bone profile was iter¬ 
atively eroded until it was only one pixel thick using a 
Hilditch skeletonisation procedure,!^ and nodes and ter¬ 
mini were automatically detected by inspecting the local 
3x3 neighbourhood. If the centre pixel of the 3x3 
neighbourhood is a skeleton pixel and one and only one 
of the 8 other pixels is a skeleton pixel the centre pixel 
is classified as a terminus (Tm) indicating that the strut 
ends in this pixel. If the centre pixel is a skeleton pixel 
and three or more of the 8 adjacent pixels are skeleton 
pixels the centre pixel is classified as a node (Nd) indi¬ 
cating that two or more struts join in this pixel. Comp- 
ston et al. have argued that the ratio between nodes and 
termini (Nd/Tm) is an expression of the connectivity of 
the trabecular network.l^ Consequently, only the node- 
to-terminus ratio (Nd/Tm) from this analysis of the his¬ 
tological sections is used in the present study. 


A. Basic definitions 


We start with the algorithm to find strands in each image. 
A strand is a connected trabecular path, i.e., a chain of 
one or more struts, with the whole path going in approxi¬ 
mately the same direction. In our algorithm, we use eight 
directions, labelled by the points of the compass: north 
(N), northeast (NE), east (E), southeast (SE), south (S), 
southwest (SW), west (W), or northwest (NW). “North” 
is the anterior direction, which corresponds to “up” in 
the images shown in this article. A node is a pixel that 
is joined to strands in at least three of these eight di¬ 
rections, any two of which are at least 90 degrees apart 
(e.g. N, E, and SW; but not N, NE, and E). The node 
strength of a pixel is 0 if it is not a node at all, but oth¬ 
erwise depends on the lengths of the strands that meet 
at the node and the pQCT values of the pixels in these 
strands. 


B. CT values of bone 


These basic definitions could be implemented in a variety 
of ways, depending for example on the definitions of “con¬ 
nected” and “same direction”, and on how the pQCT 
values of the pixels are used. In our algorithm, the first 
step is to remove marrow from further consideration, by 
applying a bone threshold filter. The thresholded pQCT 
values will be referred to in this section as CT values of 
bone, denoted b. Thus if a is the pQCT value of a pixel, 
then 


b = 


a - a^hreshold -f 

0 otherwise. 


a > a^hreshold 


(1) 


We choose a^^reshoid _ 275 ^ corresponding to a BMD 
value of 24 mg/cm^. This is the soft tissue threshold used 
in symbol encoding for complexity measures by Saparin 
et al.l^The threshold was determined as the pQCT value 
of the densest non-bone tissue tested, plus the standard 
deviation of the noise (a calibration protocol for finding 
the threshold based on a phantom will be developed in 
the future). 


C. Strands 


In order to explain the central algorithm, we begin by 
considering a fictitious image consisting of just one row. 
The CT values of the pixels in the row form a sequence 

^ 0 : bi,b2,bs, • • • . 

We will define a corresponding sequence of strand 
strengths, 


So, Si, S2, S3, • • • , 


where each Sj describes the pattern of densities to the 
left of pixel j in the row. We begin by setting sq = 0, 
and then we recursively define 




V ^-1 J 


( 2 ) 


where T is a transmission eonstant between 0 and 1. 
Note that if all of the CT values have a common value b 
then 


+ Tsj-i) — {b + T {b -^Tsj-2)) 

= b{l+T)+T^Sj-2 = --- 
= b{l + T+--- + T^), 

since sq = 0. As j increases, Sj approaches an upper 
bound of Thus, for the simple case of a uni¬ 

formly dense row, we have defined the strand strength 
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to be -—— times the CT values of bone. This constant 

-—— can be interpreted as a characteristic length of 

the method, which we can vary depending on the length 
of strand that we believe to be most important to the 
strength of the bone. For the present study, we have cho¬ 
sen T = 0.95, corresponding to a characteristic length of 
20 pixels, i.e. 4 mm. For the general case of a variable 
CT values row, the presence of the transmission constant 
T in the recursive formula ensures that the CT value of 
pixels closest to the pixel have the greatest effect on 


Sj. Finally, the factor 


min{ 6 j, 6 j_i} 

6,_i 


ensures that Sj de¬ 


pends strongly on hj^ so that any weak link in a chain of 
high-density pixels lowers the strand strength dramati¬ 
cally. 


The difficulty in generalising this definition to a 2- 
dimensional image is of course that there are infinitely 
many directions but only four that are parallel with the 
pixel grid. We have resolved this in a practical but ad 
hoc fashion. As already stated, we consider eight direc¬ 
tions. Consider first the definition of strand strength in 
the leftwards (W) direction. ^From pixel (i,j) (at row 
i and column j) we consider that there are five length-3 
paths leading approximately leftwards: 


(i) {i -h 1 , j — 2 ) ^ (i + 1 , j — 1 ) ^ 

(ii) {i + l,j -2) ^ {i,j -1) ^ 

(iii) {i,j-2) ^ {i,j-1) 

(iv) -2) ^ {i,j-1) ^ 

(v) {i -l,j -2) ^ {i- l,j - 1 ) ^ 


as illustrated in Fig. 

<i) j-2 j-1 j <ii) j-2 j-1 


<iii) j-2 7-1 7 





FIG. 2. The five possible length-3 paths in a leftwards 
(“west”) direction from pixel {i,j). The start pixel (i,j) is 
marked in black. 


D. Strand strength 


and j — 2. We then apply the recursive formula 


Sj — I bj—i + T (bj—2 + Tsj— 2 ) 


bi- 


’i -2 


(3) 


f rain{bj,bj-i} 

\ 


which is the formula in Eq. (i), applied twice. (We must 
now set Si = 5o = 0.) We now have five Sj values, 
one for each of the five paths, i.e., ..., 

We multiply all but the third (Fig. I^iii)) of these 
five values by a bending coefficient 7 between 0 and 
1 , to penalise strands that bend away from the hor¬ 
izontal direction. Then the maximum of these five 
numbers is our leftwards (W) strand strength = 

max|7s}7 75} 7 s) 7 js) 7 js) '}. 

Strand strengths in the other seven directions 
, Sjj , , • • •, are determined in a similar 

manner. For diagonal directions, the geometry of the 
five length-3 paths used in the computation is slightly 
different, so a different bending coefficient 7 is used. We 
chose 7 = 0.6 for the horizontal and vertical directions 
and 7 = 0.8 for the diagonal directions. These constants 
were chosen so as to give mean strand strengths that 
are approximately invariant to image rotations by 
an arbitrary angle, which were verified by numerical 
experiments on tibial images. 


E. Node strength 


Finally, we calculate the node strength aij at each pixel. 
There are k = 16 possible ways of choosing 3 out of 
the 8 directions in such a way that each pair of chosen 
directions makes an angle of at least 90 degrees. For 
example, E, N, and SW comprise an allowable choice, 
but E, N, and NE do not. At each pixel, for each al¬ 
lowable choice of three directions, we calculate the min¬ 
imum of the three strand strengths in these directions, 
e.g., = min{5'^^, The node strength (Tij 

is the maximum of the 16 minima subtracted by a 
minimum strength constant 

Gij = ^ max ^ } - 5o. (4) 

Pixels with a positive node strength are called nodes. The 
mean node strength over the region of interest (ROI) is 
called the node strength of the region, 

NdStr=^ E 

ijGROI 


Eor each of these paths, we define 6 ^, 6 j_i and bj -2 Rs the The purpose of subtracting the minimum strength con- 
CT values of the pixels in the path in columns j^j — 1, stant Sq is to allow us to ignore short “strands”, which 
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FIG. 3. Node strength algorithm illustrated on a region of trabecular bone taken from the full section shown in Fig. (a) 
original image; (b) horizontal strands; (c) vertical strands; (d) diagonal strands (“northwest - southeast”); (e) diagonal strands 
(“northeast - southwest”); (f) node strength. 


are often just transverse sections of trabeculae with an 
apparent width of more than one pixel. The trabeculae 
in our images have an apparent width of approximately 
two pixels, or 0.4 mm, which is higher than the true av¬ 
erage trabecular width, due to the 1 mm thickness of the 
CT image and the 0.2 mm pixel size. Thus, we wish to 
ignore strands of bone that are only two pixels long. The 
strength of such a strand depends on the pQCT values of 
the two pixels. Following Eq. the strength s of such 
a strand with CT value bo = bi = b is s = b. Based on 
examination of the resulting node strength images (such 
as Fig. (f)), we find ad hoc a minimum strength con¬ 
stant of 225 (corresponding to a CT value of 500, Eq. 0, 
corresponding to a BMD value of 331mg/cm^, which is 
higher than the mean CT value for trabecular bone). So 
strands of lower strength than this will be ignored by the 
algorithm. 

We note that our complete algorithm involves two thresh¬ 
olding steps: one at the beginning, when pQCT val¬ 
ues are converted to CT values of bone; and one at the 
end, when the minimum strength constant is subtracted. 
These are not equivalent to one thresholding step with 
a higher threshold. The purpose of the first threshold 
filtration is to ignore marrow. The purpose of the second 
threshold filtration is to ignore short strands, and it can 
be seen as an alternative to skeletonisation. 

Numerical experiments have shown that relative (cross¬ 
subject) values of NdStr are stable with respect to the 
choice of the bone threshold value (( 2 ^hreshoid^ min¬ 
imum strength constant. Specifically, if either of these 
constants is varied by ±10% from the values used in this 
article, and NdStr is recomputed for all subjects, the re¬ 
sulting NdStr values are very strongly linearly correlated 
with the original NdStr values (Pearson correlation co¬ 


efficient r > 0.995 in all cases). The absolute values of 
NdStr depend on the choice of these constants: increasing 
the bone threshold by 10% leads to a 40% mean decrease 
in NdStr (percent decrease averaged over all subjects), 
while increasing the minimum strength constant by 10% 
leads to an average decrease in NdStr of 11%. Never¬ 
theless, and as already mentioned, the relative (cross¬ 
subject) variation of NdStr remains constant. 

A key property of NdStr is that it depends on both the 
geometry of the trabecular network and the CT values 
of the trabeculae. NdStr depends linearly on thresholded 
CT values, if the CT values of all pixels are varied by the 
same factor. It follows that images with the same geom¬ 
etry but different BMD will have different NdStr values. 
On the other hand, two images with the same BMD but 
different geometry can have different NdStr values. This 
is apparent from the description of the method, and con¬ 
firmed by Eigs.[^and[^ as discussed below, and also by 
the results discussed in Section lYl 


F. Illustration 

To illustrate the analytical method, we now present the 
results in visual form for an enlarged region near the 
bottom (posterior) of the slice in Eig. The enlarged 
region of the original image is shown in Eig. ^a). Parts 
(b) to (e) of Eig. show directional strand strengths, 
and part (f) shows the final node strength plot. Each 
directional strand strength plot shows the sum of two 
strand strengths at every pixel, in opposite directions: 
east/west; north/south; northwest/southeast; and north¬ 
east /southwest. In each of the directional strand strength 
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FIG. 4. Trabecular micro-architecture of a section with BMD 107mg/cm^. (a) original image; (b) the node strength of each 
pixel. The mean node strength is 71.2. 


plots, the strands in the given direction are shown with 
the highest intensity, but most of the trabeculae are still 
visible, even if faintly. In contrast, in the node strength 
plot (part (f)), most of the trabeculae are invisible. This 
is because of the subtraction of the minimum strength 
constant. In this example, there are almost no nodes in 
the right half of the image. This correctly describes the 
micro-architecture of the original image in that region, 
which contains many trabeculae but few that cross each 
other to make a lattice-like micro-architecture. The left 
half of the image contains many nodes. Notice that, in 
the node strength plot, the nodes seem to be thicker than 
in the original image. This is because the trabeculae in 
the original image are actually slightly thicker than they 
appear, the outer pixels being dimmer (i.e. lower CT val¬ 
ues) and thus not easily registered by the eye. Since the 
outer pixels near the apparent nodes in the original image 
are almost as well-connected as pixels in the centres of 
the nodes, they have large node strengths, and are very 
visible in the node strength plot. 


The two specimens depicted in Figs. andhave compa¬ 
rable bone mineral densities, but their trabecular micro¬ 
architecture are visibly different. In Fig.[^ the specimen 
has a trabecular BMD of 107 mg/cm^, which is near the 
median value (97 mg/cm^) of the specimens in this study. 
On the left is the original image, and on the right is the 
node strength plot. Notice that there are a lot of nodes 
in most of the outer areas, with the notable exception of 
a region near the bottom left. The mean node strength 
is 71.2. 


The specimen in Fig. has a trabecular BMD of 94 
mg/cm^, which is only 12% lower than that of the speci¬ 
men shown in Fig.|^ but it has substantially fewer nodes. 
The mean node strength is only 42.2, which is 40% lower 
than that for the specimen shown in Fig. This reflects 
the lack of a strong lattice-like micro-architecture in the 
original image. 


IV. RESULTS 

The long range node strength of each pixel of each image 
was computed as described in the previous section, and a 
mean node strength, NdStr, was calculated for each im¬ 
age. We also computed the bone mineral density (BMD) 
from the pQCT slices, as explained in Section [ITU A scat¬ 
ter plot showing NdStr versus BMD for all 26 specimens 
is shown in Fig.[^ There is a strong positive correlation, 
which we quantified in three ways. Firstly, Pearson’s cor¬ 
relation coefficient is r = 0.87, indicating a very strong 
linear correlation. Secondly, since the scatter plot clearly 
suggests a nonlinear relationship, we fitted an exponen- 
tiai curve to the data and then found Pearson’s correia- 
tion coefficient to be r = 0.97. Thirdiy, the Spearman 
rank correiation is p = 0.98, indicating a very strong cor- 
reiation. The Spearman rank correiation coefficient is a 
robust nonparametric correiation measure that is appro¬ 
priate when littie is known about the distributions and 
nature of the correiation between the variabies.^^ 

We aiso compared NdStr with the node-terminus ra¬ 
tio, Nd/Tm, in the topoiogicai node-strut anaiysis intro¬ 
duced by Garrahan et ai.^ The two measures are simi- 
iar in phiiosophy, because they both quantify the nodes 
in the trabecuiar network. However, the definition of 
Nd/Tm is highly localised: after the skeletonisation pro¬ 
cess has eroded the trabecular network to a thickness of 
one pixel, each pixel is classified as node, strut, or termi¬ 
nus depending on its 3 x 3 grid of nearest neighbours (in- 
ciuding the originai pixei). For the present study, the im¬ 
ages used to compute Nd/Tm have a pixei size of 10 /im, 
so the ciassification is made on the basis of a 30 pm x 
30/im region. In contrast, node strength is semi-giobai, 
taking into account ionger strands to a degree controiied 
by the transmission strength constant. In the present 
study, with this constant set to 0.95, the method has a 
characteristic iength of 4 mm, in the sense described in 
Section im 

The node-terminus ratio, Nd/Tm, was caicuiated using 
histomorphometry performed on bone biopsies as de¬ 
scribed in Section |III[ Recall that these biopsies were 
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FIG. 5. Trabecular micro-architecture of a section with BMD is 94mg/cm^. (a) original image; (b) the node strength of each 
pixel. The mean node strength is 42.2. 


from the same regions of the same donors as the pQCT 
slices from which NdStr has been computed. Pearson’s 
correlation coefficient for the relationship between NdStr 
and Nd/Tm is r = 0.62, and the Spearman rank correla¬ 
tion coefficient is p = 0.65. We also measured the corre¬ 
lation between Nd/Tm and trabecular BMD: Pearson’s 
correlation coefficient is r = 0.64, and the Spearman rank 
correlation coefficient is p = 0.61. Using either measure 
of correlation, mean node strength is more strongly corre¬ 
lated with trabecular BMD than the node-terminus ratio 
is correlated with either of these variables. 

NdStr 


150 


100 


50 



20 40 60 80 100 120 140 BMD 

FIG. 6. Scatter plot of mean node strength vs. trabecular 
bone mineral density (mg/cm^). 


V. DISCUSSION 


We have introduced a new morphometric measure for 
characterising the micro-architecture of trabecular bone, 
long range node strength^ which measures the degree to 
which a pixel in a 2-dimensional bone image has long- 
range connectivity in three or more directions, each at 
least 90 degrees from the others. In addition, we have 
calculated the mean node strength, NdStr, for each of 
the 26 bone samples considered in the study. We have 
found that NdStr has a strong positive correlation with 


trabecular BMD (r = 0.97, after exponential transforma¬ 
tion). Furthermore, we have ascertained a strong corre¬ 
lation (r = 0.62) between NdStr and the established his- 
tomorphometric measure, node-terminus ratio (Nd/Tm). 
Moreover, qualitative comparison of images with similar 
trabecular BMD but different mean NdStr (see Figs. 
and suggest that NdStr successfully quantifies how 
“lattice-like” the micro-architecture is. 

Further studies, including either clinical data or mea¬ 
sured bone strength, are needed in order to determine 
the utility of this measure relative to other existing mor¬ 
phometric measures. Such studies could also determine 
the most useful choice of the parameters that appear 
in the algorithm: the transmission constant, minimum 
strength constant, and bending coefficients. For example, 
these constants could be chosen to maximise correlations 
with bone strength. In the future, the sensitivity of the 
method on these parameters and on different CT settings 
(e.g., CT pixel size, slice thickness, mean CT value of 
bone) should be systematically investigated. Moreover, 
a future study on synthetic data could be used to verify 
that NdStr measures structure and not just bone volume 
or mass, and specifically that it finds nodes. Also differ¬ 
ent skeletal sites should be analysed with our approach 
in order to test its potential to describe structural differ¬ 
ences. 

We would like to note that our algorithm for computing 
node strength is dependent on pixel size, and is thus un¬ 
suitable for absolute comparisons between studies. How¬ 
ever, this limitation is not unique to the node strength 
measure. For example, Guggenbuhl et al have showed 
using CT images with different thickness (1 mm, 3 mm, 
5 mm, and 8 mm) that the outcome of texture analy¬ 
sis depends substantially on the slice thicknesJ^. In the 
present study we have used pQCT equipment with an 
in-plane pixel size of 200 pm x 200 pm and a slice thick¬ 
ness of 1 mm. We have not conducted a formal inves¬ 
tigation into the influence of slice thickness on the long 
range node strength, but it is fair to assume that the 
long range node strength is similarly affected by the slice 
thickness. If further studies were to confirm the practical 







value of the measure, an implementation could be devel¬ 
oped to produce node strengths that are broadly compa¬ 
rable between images with different pixel size. However, 
the technological development of high-resolution pQCT 
equipment is progressing very quickly, and already high- 
resolution pQCT scanners exist that can image a tibia at 
an isotropic pixel size of approximately 100 /im, which 
is comparable to the trabecular thickness in the human 
proximal tibisP^l, thus making the slice thickness less of 
an issue. 

Some previous studies have investigated the trabecular 
bone micro-architecture using texture analysis applied 
on X-ray images of bone^^HMl^ However, these tech¬ 
niques are not comparable to the method presented in 
the present study as they are based on projections of 
a 3-dimensional trabecular network on a 2-dimensional 
plane, whereas our method is applied to 2-dimensional 
sections obtained through the 3-dimensional trabecular 
network. Nevertheless, Ranjanomennahary et al. very 
recently showed some significant correlation did exist be¬ 
tween radiograph based texture analysis and fiCT based 
unbiased 3-dimensional measures of trabecular micro- 
architectur^^. Apostol et al compared 3-dimensional 
measures of micro-architecture based on 3-dimensional 
synchrotron radiation CT with more than 350 texture 
parameters obtained from simulated radiographs that 
were created from the synchrotron radiation CT data 
setJ^. They found using multiple regression analysis 
that a combination of a subset of texture analysis param¬ 
eters correlated to the 3-dimensional measures of micro¬ 
architecture. However, they had to use at least three tex¬ 
ture analysis parameters in the correlation with each of 
the 3-dimensional measures of micro-architecture. Cortet 
et al. investiga ted C T images of the distal radius using 
texture analysi^^^^. The CT image used by Cortet et al 
used similar pixel size to those used in the present study. 
They analysed the CT images using the traditional node¬ 
strut analysis and te xture analysis including the grey 
level run length methocP^^^. As illustrated in the present 
work there is a moderate correlation (r = 0.64) between 
the node-strut analysis and the NdStr measure, which 
we believe illustrates that the NdStr measure captures 
somewhat different information from that obtained with 
the node-strut analysis. Grey level run length is based 
on runs that have exactly the same grey level, and is 
therefore very sensitive to the choice of discretisation of 
the grey levels. In contrast, the NdStr measure treats 
the grey level as a continuous parameter, and thus it is 
able to detect runs with variable grey levels, and is con¬ 
sequently less sensitive to discretisation choices. 

Another advantage of the method is that the proposed 
method utilises all of the CT value information in the 
pQCT images, as no binarization of the images is per¬ 
formed. 

In the present study we have not applied the developed 
method to the histological sections directly. The rea¬ 
son for this is that the histological sections only cover a 


limited area of interest and thereby a limited trabecular 
length, which renders the NdStr less meaningful. How¬ 
ever, this is a limitation of the histological examination 
procedure where it is only feasible to investigate a smaller 
sample (biopsy) of a larger structure (the proximal tibia) 
and not a limitation of the proposed method. In ad¬ 
dition, the histological images are segmented into bone 
and marrow by their nature and it is thus not possible 
to assign an intensity value to pixels, analogous to a CT 
value, which is needed by the algorithm in its present 
form. However, it is possible to apply the concept of 
long range node-strut analysis to such binary images but 
this is outside the scope of the present investigation. 

In the present study the long range node-strut analy¬ 
sis was applied to pQCT images of the proximal tibia. 
However, we would like to stress that the method is not 
limited to this skeletal site and thus can be applied to 

2- dimensional CT images obtained from any skeletal site 
like e.g. the vertebral body or the calcaneous. 

Finally, we note that the general method of long range 
node-strut analysis provides more than just the mean 
node strength. In the present study, we have focused 
on mean node strength for simplicity, but the intermedi¬ 
ate measures of direetional strand strength^ used in the 
computation of node strength, may be useful in them¬ 
selves as a measure of directional strength. In the present 
study the long range node-strut analysis has been applied 
to 2-dimensional pQCT images obtained in the horizon¬ 
tal plane. However, the trabecular micro-architecture of 
the proximal tibia is mostly isotropic in the horizontal 
plane, whereas the micro-architecture in vertical direc¬ 
tion is highly anisotropic compared with the horizon¬ 
tal plane^. As jaCT scanners become more and more 
prevalent and as the pixel size and imaging capacity of 
pQCT scanners steadily improve it will be an important 
future task to generalise the long range node-strut anal¬ 
ysis into three dimensions and to apply the technique to 

3- dimensional data sets obtained from such equipment. 
In particular, the directional strand strength could be 
used to investigate anisotropic differences of the trabec¬ 
ular micro-architecture of such 3-dimensional data sets. 
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